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ABSTRACT 

An electron-positron collider operating at a center-of-niass energy Ecm can col¬ 
lect events at all lower energies through initial-state radiation (ISR or radiative 
return). We explore the capabilities for radiative return studies by a proposed 
high-luminosity collider at Ecm = 250 or 90 GeV, to hll in gaps left by lower- 
energy colliders such as PEP, PETRA, TRISTAN, and LEP. These capabilities 
are compared with those of the lower-energy e’''e“ colliders as well as hadron 
colliders such as the Tevatron and the CERN Large Hadron Collider (LHC). 

Some examples of accessible questions in dark photon searches and heavy flavor 
spectroscopy are given. 

PACS codes: 13.66.Be, 13.66.De, 13.66.Hk 

I Introduction 

An electron-positron collider operating at a center-of-mass energy Ecm can collect events 
at all lower energies through initial-state radiation (ISR). This radiative return process has 
been used to good advantage in e'''e“ colliders such as DAd>NE, PEP-H, KEK-B, and LEP 
[1-4]. In the present paper we explore the capabilities of a higher-energy high-luminosity 
collider such as that envisioned by CERN (FCC-ee) [5] or China (CEPC) [6], operating 
at Ecm — 250 or 90 GeV (functioning as a Giga- or Tera-Z factory at the latter energy) [7], 
to perform radiative return studies of physics at lower center-of-mass energies. 

In order to fairly assess the capabilities of future colliders with past and present colliders 
it is necessary to specify the total integrated luminosity expected to be collected by future 
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Table I: Projected luminosities for the CEPC [8] and FCC-ee [9]. These values are used 
throughout the text. 



= 90 GeV 

= 250 GeV 

GEPG 

0.5 ab-i 

5 ab“^ 

FGG-ee 

50 ab-^ 

10 ab“^ 


colliders. Based on current design reports, over 2 interaction points the CEPC is expected 
to collect 500 fb“^ on the Z pole, which corresponds to approximately 1 x 10^° Z’s, and 5 
ab“^ at Ecm — 250 GeV [8]. The FCC-ee, over 4 interaction points, is expected to collect 50 
ab“^, which is roughly 1 x 10^^ Z’s, at Ecm — 90 GeV and 10 ab“^ at Ecm — 250 GeV [9]. 
Table I summarizes these numbers. For convenience where the exact number of events is 
not imporant, we shall occasionally quote results for a nominal integrated luminosity of 1 
ab“^. These values may be rescaled appropriately. 

We briefly review some previous uses of radiative return in Section II. In Section III 
we study narrow resonance production, while the discussion is extended to continuum 
production in Section IV. Section V compares the reach of and hadron colliders for 
two benchmark processes: “dark photon” and bb production. Some processes of interest 
in heavy flavor spectroscopy are noted in Section VI. We conclude in Section VII. Some 
calculational checks are contained in two Appendices: VIII and IX. 

II Some previous uses of radiative return 

Considerable use has been made of radiative return in previous experiments using electron- 
positron colliders. In Table II we summarize some parameters of experiments at these 
colliders [11,12]. Maximum instantaneous luminosities of circular e’''e“ colliders are plotted 
versus year in Fig. 1. 

A KLOE at DA$NE 

The DA<hNE accelerator at Frascati operates near or at the CM energy (1020 MeV) of the 
0 resonance. It has studied the cross section —?■ 7r’''7r“ at lower CM energies via the 

process e+e“ —?■ yTr+Tr”, where the photon is emitted in initial-state radiation, with the 
main purpose of reducing the error in the hadronic vacuum polarization contribution to 
the anomalous magnetic moment of the muon. Three sets of data are reported: 141.4 
pb“^ studying the interval 0.35 < < 0.95 GeV^ [14], 240 pb“^ studying the same 

interval [15,16], and 230 pb“^ studying 0.1 < < 0.85 GeV^ [17]. KLOE also has 

searched for “dark photons” below 1 GeV decaying to and as noted in more 

detail later [18,19]. 

B CLEO at CESR 

The GLEO Gollaboration used initial-state radiation to search for the state V(3872) [20] 
in collisions at the Gornell Electron Storage Ring (GESR). The absence of a signal 
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Table II: Instantaneous and/or integrated luminosities achieved at some e'''e“ colliders. 
Based in part on Section 30 of Ref. [11], with values from Ref. [12] for PETRA, PEP, and 
TRISTAN. We thank G. Alexander and S. L. Wu for help with some of these estimates. 


Collider 

Detector 

CM energy 
(GeV) 

Max. C 
(10^° cm“^s“^) 

/ Cdt 
(fb-i) 

DA$NE 

KLOE 

1.02 

453 

2.5 



1.00 

453 

0.23 

CESR 

CLEO 

9.46-11.30 

1280 at 10.6 GeV 

15.1 

PEP-H 

BaBar 

10.58 

12069 

424.7 



10.18 


43.9 

KEK-B 

Belle 

9.46-10.89 

21083 

980 

PEP 


29 

60 

1.167“ 

PETRA 


46.8'^ 

24 at 35 GeV 

0.817“ 

TRISTAN 


64^ 

40 

0.942^^ 

LEP 


Mz 

24 

0.808“ 



> 130 

34-90 

2.980“ 


“ Summed over detectors DELCO, HRS, MAC, Mark II, TPC/27 
^ Maximum value 

Summed over detectors CELLO, JADE, Mark J, PLUTO, TASSO 
^ Summed over detectors AMY, TOPAZ, VENUS 
® Summed over detectors ALEPH, DELPHI, L3, OPAL 


served as partial evidence that the state did not have spin 1 and negative parity and 
charge-conjugation eigenvalue. 

C BaBar at PEP-II 

The initial-state radiation process has been used to great advantage by the BaBar Collab¬ 
oration at PEP-H. Just in the past three years, papers have appeared on the production of 
7r'''7r“7r“''7r“ [21]; [22]; 7r’''7r“ [23]; [24]; 

pp [25]; K^K~ [26]; [27]; and a variety of hnal states with two neutral 

kaons [28]. The CM energies and integrated luminosities in Table H are those quoted 
for BaBar in the last paper. The hnal states involving light hadrons contribute to reducing 
the uncertainty on the hadronic vacuum polarization contribution to and to the running 
of the hne structure constant in precision electroweak studies, while those involving J/'i/ 
and t/(2S') are of interest for resonant structures. 

D Belle at KEK-B 

Since 2011 the Belle Collaboration has produced a couple of initial-state radiation studies, 
involving production of J/'0A'+A'“ and J/piKsKs [29]; and [30]. The focus of 

this work has been the search for resonant substructures in the hnal states. 
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Figure 1: Maximum instantaneous luminosities of circular e'*'e colliders versus time. 
Adapted from Fig. 1 of Ref. [13]. 

E LEP Collaborations 

The LEP entries in Table II refer to total integrated luminosities at various energies. Results 
from specific detectors are not always based on these totals, as many of them were obtained 
before the full data sample was available. 

E.l ALEPH 

The reaction 7 isr has been used by the ALEPH Collaboration [31] to study 

the cross section and forward-backward asymmetry for —)■ in the CM energy 

range 20-136 GeV. 

E.2 DELPHI 

The helicity structure in e’''e“ —)■ ii~ is particularly sensitive to new physics in the CM 
energy range around 80 GeV, which is accessible through the process e+e” The 

DELPHI Collaboration [32] has studied this process. Ending no evidence for new physics. 

E.3 L3 

The L3 Collaboration has used ISR to measure muon pair production in collisions 

between 50 and 86 GeV [33]. This study was motivated in part by the need to fill a gap 
between the maximum energy of the TRISTAN accelerator at KEK (about 62 GeV) and 
the mass of the Z. 
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E.4 OPAL 


The OPAL Collaboration pioneered the use of ISR at LEP to £11 the aforementioned gap 
between 62 GeV and Mz [34]. No deviations from the standard model were found, albeit 
with a very early data sample. The collaboration used radiative fermion pair events to 
perform a LEP beam energy measurement [35]. Fig. 1 of this reference gives an idea of the 
yield that may be expected from radiative return studies with beam energies approaching 
200 GeV. An extensive study was made, accumulating a total of 1.132 fb“^ at LEP CM 
energies between 183 and 207 GeV. Aside from a copious Z peak, the numbers of events 
per 2 GeV subenergy at a subenergy of 125 GeV were about 200 for qq, a dozen for 
and three dozen for r+r“. Supposing one had a sample of 1 ab“^ at Eqm = 250 GeV, one 
might expect 0{2 x 10^) hadronic events, 0(10^) events, and a few tens of thousand 

T~^T~ events, per 2 GeV bin in invariant mass around 125 GeV. 


Ill Resonance production 


The cross section for electron-positron production of a vector meson resonance R with mass 
mji and partial width Pee decaying to a final state / with partial width P/ may be 
written near resonance as 


cr(e^e —)■ R —)■ /; s) 


127rPeeP f 

{s - m\y + (mijPij)2 


( 1 ) 


where s = and and P/j are the resonance mass and total width. ^ For the T(4S'), 

whose decays are almost exclusively to BB final states, the leptonic branching ratio is 
quoted by the Particle Data Group [11] as 1.57 x 10“^ while the total width is 20.5 MeV, 
leading to a leptonic partial width Pee = 0.322 keV. We shall use this value, noting that 
it is mildly inconsistent with the Particle Data Group’s average of 0.272 keV. The mass 
is 10.5794 ± 0.0012 GeV; the cross section at the resonance peak is about 2.06 nb. The 
resonance shape is shown at the left in Fig. 2. 

A resonance R may be produced by the radiative return process e’''e“ —?■ yi?, where the 
electron or positron of beam energy E = Ecm/‘^ radiates a fraction 1 — x of its energy and 
is left with energy xE. Neglecting the small electron mass, the squared effective mass of 
the e^e“ system is then xs. An electron beam of energy E radiates a photon and ends up 
with an energy xE with a probability per unit x [37] denoted by 


/e(x, ^/s, Pt, cut) 


a 1 4- , E 

-In- 

TT 1 - X Pt, cut 


( 2 ) 


where the minimum photon transverse momentum Pt, cut provides a collinear cutoff.** In 
the absence of an explicit choice of cutoff, it is provided by the electron mass rrie, which we 
shall use in much of what follows. The cross section for production of the resonance R by 
radiative return, where R decays to the final state /, is then 


cr(e’'"e ^ yR —)■ 7 /) = — In— [ dx ^ cr(e'*~e —)■ i? ^/; xs) , (3) 

TT me Jo 1 — X 

II An extensive discussion of possible modifications of this expression, including multiplication of F by 
the factor sjrrdp to ensure 1/s behavior of the cross section at high s, is given in Ref. [36]. 

**The numerator of the logarithm is sometimes taken to be 2A = %/s,. 
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Figure 2: Cross section for e+e —)■ T(4S') (left) and including the emission of a photon at 
an e^e“ collider with CM energy Eqm (right). 

where the factor of two comes from the fact that either lepton can radiate the photon. In 
the narrow-resonance approximation, the integral in this expression can be done in closed 
form, with the result 

cr(e^e“ —)• yi? 7/) ~ 24«7rln-——^ ^ , (4) 

rrie ^ — Xq rriR s 

where Xq = s and Bj = Tj/Tr denotes the branching fraction into the hnal state 
/. The cross section for e’''e“ —)■ T(4S') including the emission of a photon is shown as a 
function of CM energy in Fig. 2 (right). 

The proposed high-energy electron-positron colliders at CERN and in China anticipate 
integrated luminosities of 50 ab“^ and 0.5 ab“^, respectively, at CM energy 90 GeV, and 
10 ab“^ and 5 ab“^, respectively, at 250 GeV [8,9]. The observation of a new resonance 
with at least 10 events would then require cross sections of at least 0.2 and 20 ab at CERN 
or China, respectively, at 90 GeV, or at least 1 and 2 ab, respectively, at 250 GeV. 

Fig. 3 illustrates contours of equal cross section for an e’''e“ collider with CM energy 
90 (top) and 250 (bottom) GeV to produce a resonance of mass rriR via radiative return. 
These results imply a cross section of 9.17 fb for the T(4S') produced by radiative return 
at Ecm = 90 GeV, given an assumed leptonic partial width of Fee = 0.322 keV [11]. For a 
given Ecm, the lowest sensitivity appears to occur for a resonance mass roughly equal to 
Ecm/‘2^, he., the beam energy. 

The results of Fig. 3 can be expressed in more universal form. In the narrow-resonance 
approximation the predicted radiative return cross section, Eq. (4), is directly proportional 
to TeeBf, so the ratio cr(e’''e“ — ?■ jR —)■ 7/)/Feei3/ is a function only of s and itir. In Fig. 4 
we plot this ratio as a function of resonance mass for two values of Ecm- 

The leptonic widths of known and hctitious quarkonium states can serve as benchmarks 
for the interpretation of Figs. 3 and 4. They are summarized for IS states in Table III 
[11,38]. (The bound state of an actual top quark t of mass ~ 173 GeV/c^ and a 7 is highly 
unstable due to the weak decay of the t or t.) 
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Figure 3: Contours of equal cross section for radiative return production of a resonance with 
leptonic width Fee (assuming 100% branching fraction to a hnal state /). Top: Ecm = 90 
GeV; bottom: Ecm = 250 GeV. 


Table III: Leptonic widths of known and hctitious IS quarkonium states. 


IS 

state 

Quark 

charge 

Quark 
mass (GeV) 

Tee 

(keV) 

Ref. 

Jl^ 

2/3 

1.4 

5.55 ±0.14 ±0.02 

[11] 

T 

”1/3 

4.8 

1.340 ±0.018 

[11] 

Toponium 

2/3 

40 

6.5 

[38] 

(hctitious) 

2/3 

45 

6.7 

[38] 
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The leptonic width of an S-wave quarkonium QQ bound state of mass M is given by [39] 

ree(QQ) = ^^^l^(0)P, (5) 

where cq is the quark charge, |\I/( 0 )p is the square of the wave function at the origin, and 
we have neglected relativistic and QCD corrections. One sees that the leptonic widths in 
Table III are governed primarily by the square of the corresponding quark charge. Over a 
wide range of quark mass, the decrease of the 1 /M^ factor is approximately compensated 
by a corresponding growth of |\l/( 0 )p. 

Note that for weakly coupled resonances typically the width is oc M, which needs to 
be accounted for when interpreting Figs. 3 and 4. In this case the factor of Tee/M in 
Eq. (4) is mostly independent of mass and the cross section is governed purely by the 
factor (1 + Xq )/(1 — Xq), falling monotonically with decreasing resonance mass. 

This discussion is oversimplihed because it neglects the off-shell process e’''e“ jZ* —)■ 
'jR, in the absence of assumptions about how R couples to Z. However, it gives an idea of 
the orders of magnitude necessary to hnd a previously missed resonance in the mass range 
below that in which the Z* contributes appreciably (e.g., below about 60 GeV, the CM 
energy accessible to TRISTAN). 


IV Continuum production 


An important quantity is the effective luminosity of a high-energy collider for studying 
any given process at lower center-of-mass energy. Dehning a{s) = a{e'^e~ 7 /; s) and 

a{s) = cr(e+e“ /; s) for a given hnal state /, the relation between the two is 


da{s) 

dx 


2al + x"^ 
71 1 — X 


In —a{s) 
Tile 


1 


( 6 ) 


where x = s/s. The subsystem CM energy may be denoted Ecm = Vs- The cross section 
per unit Ecm times an interval A of Ecm is then 


da{s) 

dEcM 


A 


AaEcM 1 + 

71S 1 — X 

Lfd{s) , 


E 

A In — (j(s) 

nip 


(7) 


where Lf is the fractional luminosity per Ecm bin of size A. Examples of this function for 
a bin width of A = 1 GeV are shown in the top curves of Fig. 5. 

For low Ecm one may take (1 -|- x "^)/(I — x) ~ 1 in Eq. (7). Integrating from Effff = 10 
GeV to E'ff^ = 30 GeV, one then Ends 

L, = In T (8) 

1 1 S TTlj^ 

For Ecm = (90,250) GeV we find Lf = (5.22,0.74) x 10“^. For a total of 1 ab“^ at 
Ecm = (90, 250) GeV this then provides a total integrated luminosity of (5220,740) pb“^ 
in the range 10 < Ecm A 30 GeV. This exceeds integrated luminosities at PEP or PETRA 
(see Table II). 
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Figure 4: cr(e~^e —)■ jR —>■ 'jf)/TeeBf as a function of resonance mass for Ecm = 90 (left) 
and 250 (right) GeV. 




Subsystem energy (GeV) 


Figure 5: Fractional luminosity Lj as a function of subsystem energy Ecm for Ecm = 90 
(left) and 250 (right) GeV. Top curves: No minimum angle; infrared cutoff provided by 
ln(F'/me) [Eq. (7)]. Lower curves, top to bottom: Oq = 10, 20, 30, 40°. 











Integrated Luminosity 



Figure 6 : Integrated luminosity from past low energy e’''e“ colliders at their nominal center- 
of-mass energies compared to the effective luminosity through radiative return from future 
e’''e“ colliders at y/s = 90 or 250 GeV (no minimum angle; see Fig. 5 for effects of minimum 
angles). The FCC-ee curves assume an integrated luminosity of 50 ab“^ at 90 GeV and 
10 ab“^ at 250 GeV. The GEPG curves assume an integrated luminosity of 0.5 ab“^ at 90 
GeV and 5 ab“^ at 250 GeV. Integrated luminosities of PEP-II and Belle (Table II) exceed 
those achievable by radiative return at FGG-ee or GEPG running at 90 or 250 GeV. 


Given the concept of fractional luminosity we can compute the effective luminosity 
gathered at each center-of-mass energy via radiative return. This is shown in Fig. 6 for 
y/s = 90 and 250 GeV compared to the luminosity collected directly at various other 
colliders. In this hgure it is illustrated clearly that a high-luminosity high-energy e^e“ 
collider both competes with and hlls in gaps left by previous colliders. 

Gleaner signals for radiative return may be obtained at the expense of recorded events 
by demanding that the radiated photon make a minimum angle 6o with respect to the beam 
axis. Let 9 be the polar angle of the radiated photon, and 2 ; = cos 6 ', zq = cos 6 ^ 0 . Using 
Eq. (8) of [40], we find the angular distribution of the radiated photon for rUe = 0 is given 
by 


dEcMdz 


AaEcM 

TIS 


1 — z'^ 


— X 

_ 4 


(1 + ^^) + 


X 

1 — X. 


(9) 


This may be integrated between the desired limits of 6 ^, with the result 



d^a{s) 

dE(j]^dz 


AaEcM 


-cr s 


TTS 




1 + 2^0 
1- Zo 



X 1 + Zo' 

-In- 

1 — X 1 — Zq_ 


( 10 ) 


The ratio between the left-hand side and (j(s) is again a fractional luminosity and is shown 
by the lower curves in Fig. 5, again for a bin width of 1 GeV. In the limit of small 9q zz 
Pt/ the leading-logarithmic term of Eq. (10) reduces to the form in Eq. (7). 

Note that our computation of fractional luminosity utilizes factorization in the collinear 
limit. In the Appendix, we perform the exact calculation for the process e’''e“ —?■ pi^ii~ 
and hnd good agreement with our results in Fig. 5. 
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V Sensitivities of colliders for benchmark processes 

We estimate the reach of radiative return studies using electron-positron colliders for two 
benchmark processes: “dark photon” searches and bb production. We compare these sensi¬ 
tivities with those of hadron colliders. Although the latter have an advantage in total rate, 
it can only be realized with considerable background suppression, such as provided for bb 
production by the VErtex LOcator (VELO) in the LHCb experiment. 


A Dark photon search 

In this section we compute the reach for dark photons using radiative return as a concrete 
example of a search for weakly-coupled resonances. This has been previously computed 
for GeV-scale dark photons using low energies colliders like PEP-II in [41]. In this work 
we focus on the lO’s to lOO’s of GeV scale, as discussed in [42]. Other relevant work 
includes [43-51]. 

For simplicity, we assume that a “dark photon,” denoted by Z', is kinetically mixed 
with a hypercharge gauge boson with amplitude e 


r — _1 _ 1 I , 1 

A 2c, 


/2 


+ -MlZ^ 


( 11 ) 


where c„ is the cosine of the Weinberg angle and the hats denote states that are not mass 
eigenstates. After diagonalization one hnds a single massless state, identihed to be the 
photon. The would-be standard model Z and dark photon Z' also mix due to electroweak 
symmetry breaking. The mixing formulas can be worked out analytically but are not shown 
here (see [42] for the full expressions). 

The dark photon inherits couplings to fermions both from mixing with hypercharge 
and mixing with the Z. In the limit e <C 1 and Mz' -C Mz the dark photon couplings to 
fermions become photon-like and the partial width simplihes to 

T(Z' ff) = {BM'j _ (12) 


where there are Nf. colors of / with charge Qf and mass m^, and 



Am'j 


(13) 


Ignoring all quark masses except and assuming the top-antitop channel is closed, the 
branching ratio of Z' into (a convenient and low-background hnal state) is 


B{Z' p+/i-) 


3 



A(3 - 
2 


-1 


(14) 


When Mz' ~ Mz the dark photon couplings become Z-like and when Mz' ^ Mz they 
become 5-hke. This can be seen in Fig. 7 where we show the branching ratios, assuming 
the dark photon decays entirely into standard model particles. These are computed using 
e = 5 X 10“^ although for e -C 1 the branching ratios are independent of e. For simplicity, 
we only use the perturbative calculation. For low Z' masses, i.e., below a few GeV, it is 
necessary to consider threshold effects, QGD corrections, and hadronic resonances. 


11 








Z' mass (GeV) 


Figure 7: Dark photon branching ratios. These are computed using e = 5 x 10 ^ although 
for e 1 the branching ratios are independent of e. 

A.l Leptonic production 

One recent use of radiative return has been the search by the KLOE Collaboration at the 
DA<I>NE e+e“ collider for a “dark photon” U decaying to /r■*■/!“ [18] or [19], produced 
by the reaction 7 isr^ at an initial CM energy of about 1 GeV. We make projections 

for a similar analysis that can be performed for initial CM energies of 90 and 250 GeV. 

To compute the reach we consider the background to be where the 

muons come from an intermediate 7 * or Z. The search then proceeds by counting the 
number of events in the dimuon invariant mass spectrum. Since the dark photon width is 
very narrow the best significance is achieved by binning as narrowly as possible around the 
targeted Z' mass. The smallest invariant mass bin is determined entirely by the detector’s 
resolution. Typically detector resolution for muon-based searches gets worse at higher 
momentum (and equivalently higher invariant dimuon mass). We take the mass resolution 
to be given Am = m^/( 10 ^ GeV) by estimating based on the specihcation A{1/pt) = 
2 X 10“^ GeV“^ outlined in future detector designs [52]. 1 For reference, this equates to 
Am = 100 MeV for m = 100 GeV. The limit on e scales as so a 4 times increase 

in resolution only results in a 40% increase in reach on e. 

The results are shown in Fig. 8 . The current and projected limits from electroweak 
precision data (EWPT) were computed in [42]. Due to some mild tension in the electroweak 
fit in present data [53] the inclusion of a dark photon with Mz' > Mz actually improves the 
fit, which is the reason that projected precision electroweak limits are weaker. Alternative 
projections that assume the electroweak precision data converges are presented in [42]. At 

^See Appendix A in [41] for an estimation of BaBar’s mass resolution. They find the mass resolution 
to grow quadratically with mass. 
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Figure 8 : Dark photon limits at 95% C.L. on the hypercharge mixing e as a function of dark 
photon mass. The i/i = 90 GeV and 250 GeV lines show our projections with future e+e“ 
colliders with integrated luminosities specihed in Table I. Electroweak precision constraints 
(EWPT) and direct searches are taken from [42], The 100 TeV projection assumes an 
integrated luminosity of 3000 fb“^. 

masses below the current direct searches are originally taken from [54] which uses the 
Drell-Yan process pp ^ Z' ^ normalized to 7 TeV LHG data [55] to compute limits 
using the full 7 and 8 TeV LHG data set.^ For direct searches for masses above the 
limits are originally taken from [56], which recast an ATLAS dilepton search [57].^. The 
100 TeV direct searches are taken from [42] which rescaled the previously mentioned direct 
limits to 100 TeV with 3000 fb-P 

A.2 Hadronic production 

Direct photon production by hadronic collisions in the standard model proceeds through the 
subprocess qq —?■ 7 *. Assuming a “dark photon” is produced by this same process, where 
the 7 is now virtual and mixes kinetically with the dark photon Z', one can utilize Drell- 
Yan production of a lepton pair or to evaluate the sensitivity of dark photon 

tPor current direct limits on dark photons both above and below Mz we take limits from [42] rather 
than from the original studies [54,56]. 

§As far as we can tell, this limit does not stop at Mz' ~ 175 GeV for any fundamental reason, but 
rather because that is the lowest mass shown in the ATLAS results [57]. 
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Figure 9: Cross sections for Z' production [54] based on kinetic mixing at a pp collider 
with CM energy 7 TeV (red), 8 TeV (blue), and 14 TeV (green). Plotted values are to be 
multiplied by [e/(2 x lO-^)]^ = (SOe)^. 


searches in hadronic collisions. A sample calculation has been performed in Ref. [54] for 
various LHC energies; the cross sections are shown in Fig. 9. (The reach of a future 100 
TeV pp collider has been investigated in Ref. [42].) 

The peaking of parton distributions at low Feynman x favors low Z' masses. For 
example, at 14 TeV a 15 GeV Z' has a cross section of about 300 x (2500e^) pb = 750e^ 
nb. Assuming an integrated luminosity of 1 ab“^, this gives rise to at least 10 events when 
> 1.3 X 10“^^. At such a low Z' mass, however, background considerations probably 
dominate any realistic estimate of sensitivity. 

B Production of bb 

B.l Leptonic production 

The asymmetric B factories PEP-II and KEK-B have explored bb production up to CM 
energies of about 11 GeV with compelling statistics, and the upgraded KEK-B with the 
Belle-II detector will extend samples to dozens of events per attobarn. However, from about 
11 to 90 GeV the e+e“ territory is much more sparsely populated with data, as one can 
see from Table H and Fig. 6 . Radiative return studies from a Giga-Z or Tera-Z factory 
can help to hll this gap. A sample process is e’''e“ {'j*,Z*) —)■ bb, compared for direct 

production with the radiative return process —)■ ( 7 *, Z*) 71 SR —)■ bbjisR. 

We use lowest-order MadGraph [58] for our estimates of direct and radiative-return 
cross sections. For simplicity we assume that the whole PEP sample of 1.167 events per fb 
is accumulated at y/s = 29 GeV, where MadGraph predicts a(e~^e~ —)■ bb) = 36.8 pb, giving 
a total sample of about 43k events. The corresponding cross sections at 35 and 60 GeV, 
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Table IV: Comparison of direct and radiative-return e+e production of bb. 


Gollider 

Direct 
cr / Cdt 

(pb) (pb-^) 

Events 

(10^) 

E range 
(GeV) 

a (pb) 

Radiative return 

= 90 GeV = 

Evts.(10®)“ a (pb) 

250 GeV 
Evts.(103)^ 

PEP 

36.8 

1167 

42.9 

10-35 

0.494 

(24.5,0.245) 

0.066 

(660,330) 

PETRA 

25.8 

817 

21.1 

35-60 

0.410 

(20.5,0.205) 

0.039 

(391,196) 

TRISTAN 

16.2 

942 

15.3 

60-85 

12.94 

(647,6.47) 

0.256 

(2562,1281) 


“ Assuming / Cdt = (50,0.5) ab ^ at (FCC-ee, CEPC). 

^ Assuming / Cdt = (10,5) ab“^ at (FCC-ee, CEPC). 

relevant for PETRA and TRISTAN, are 25.8 and 16.2 pb, respectively. With integrated 
luminosities of 817 and 942 events per pb (Table II), one then has respective samples of 
21.1k and 15.3k events from the direct process at PETRA and TRISTAN. 

For radiative return we consider samples integrated over the Eqm ranges [10,35], [35,60], 
and [60,85] GeV, applying Eq. (7) and recalling that the beam energy E is \/s/2.' The 
results are compared with the direct process in Table IV. 

Even if all the data from PEP, PETRA, and TRISTAN are pooled, they are less than 
the sample that would be gained in the Eqm range from 10 to 60 GeV by studying radiative 
return from an e’''e“ collider at ^/s = 90 or 250 GeV. 

B.2 Hadronic production 

It is not straightforward to compare leptonic and hadronic bb production because the back¬ 
ground circumstances are different. However, the LHGb Gollaboration has demonstrated 
great sensitivity to specihc hnal states in which backgrounds can be largely overcome. As 
one example, a recent study of b hadron lifetimes [59] based on a data sample of 1 fb“^ at 
^/s = 7 TeV accumulates a sample of 229,439 ± 503 B~^ —)■ J events, with J/i^ de¬ 
caying to ii~. Given the branching fractions B{B'^ J/ipK'^) = (1.027 ± 0.031) x 10“^ 
and —)■ /i''’/i“) = (5.93 ± 0.06)% [11], this corresponds to about (3.77 x 10®)/e/ B^ 

produced, where e/ < 1 is the acceptance for the hnal state / in question. In fact, the 
production cross sections for B mesons at the LHG have been measured [60]: 


a{pp B'^ + X) 

= (38.9 ± 0.3 ± 2.5 ± 1.3)/ib , 

(15) 

a{pp B^ + X) 

= (38.1 ±0.6 ±3.7 ±4.7)/ib , 

(16) 

a{pp Bs + X) 

= (10.5 ± 0.2 ± 0.8 ± 1.0)/ib , 

(17) 


where the errors are statistical, systematic, and normalization (based on prior branching 
fraction measurements). This would correspond to about 3.9 x 10^° B^ in a sample of 1 
fb“^, yielding an estimate of e/ ~ 10%. 

^In MadGraph the logarithm is taken to be ln(Vs/PT,cut); we apply corresponding corrections of 0.9426 
and 0.9471 to the MadGraph results at ^/s = 90 and 250 GeV. 
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VI Some accessible questions in heavy flavor spec¬ 
troscopy 

Much progress in heavy flavor spectroscopy has been made using the B factories PEP-II 
and KEK-B. However, these machines were limited, as will be the KEK-B upgrade, to CM 
energies not much above 11 GeV. There are a number of questions in the spectroscopy of 
hadrons containing heavy (charm and bottom) quarks that could benefit from higher CM 
energies. Can an e'''e“ collider with CM energy 250 GeV and luminosity 10^'^ cm“^s“^ 
provide integrated luminosity to study such states significantly above what has already 
been provided by PEP, PETRA, TRISTAN, and LEP? A sharper answer can be provided 
by considering specific processes. 

A Bottomonium analogues of charmonium X, V, and Z states 

There are a number of charmonium states that appear to contain extra light quarks or to 
be admixtures of cc and charmed meson pairs. Some analogues of these have been seen 
in the bottomonium sector, but so far the X, the analogue of the V(3872), has eluded 
clearcut detection. There is an intriguing possibility that might have already been 
observed, but identified as Xfei(3P) which is close in mass and has the same quantum 
numbers [61]. Electron-positron collisions with CM energy greater than 11 GeV may be 
helpful in resolving this issue and allowing for unambiguous identification of and related 
states. 

B Pair production of narrow Bgj states 

The reaction e’''e“ —?■ Bgj + X may be used to look for the Aquark analogue of the very 
narrow D^j states seen by BaBar, GLEO and Belle [62-64]. The relevant thresholds are 
discussed in Subsection D below. 

C Doubly heavy flavor production 

With sufficient GM energy one may study such processes as 

e''“e“ — hbcc -|- X , 

e''“e“ —)■ bbbb + X , (18) 

as a precondition for producing doubly heavy mesons such as B^, B*, and doubly heavy 
baryons such as = bcq, and = bbq, where g is a light quark. Until now the latter 
have never been clearly observed, even though it is clear they must exist. As shown in 
Ref. [65], one must be able to see the (known) B^. state if one expects to be able to detect 
Sftc, so we shall estimate Be production by radiative return. We shall consider the case 
Ecm = 90 GeV, assuming that a circular e^e“ collider will spend some time as a Giga- (or 
Tera-) Z factory. 

The mass of Be is by now very well known [11]: M{Be) = 6275.6±1.1 MeV. For optimal 
Be production, one probably needs to be above B*'^Be~ threshold, which according to 
the estimate of Ref. [65] lies between 12.69 and 12.72 GeV. The cross section a{e'^e~ —)■ 
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B*~^B*~) probably rises sharply near threshold, in the same manner as a{e'^e~ —)■ 
and may be estimated as follows. 

The cross section for production of a bb pair, far enough above the T(4S'), is 

expected to be 


a{e'^e bb; s) 


47r«^ 1 

3s ' 3 


47ra^ 1 1 

3s x 3 ’ 


(19) 


and is about 180 pb at = 12.72 GeV. By comparing cross sections for and Be 
production at LHCb, Ref. [65] found the probability of a & quark fragmenting to B~ = be 
to be about 10“^. Thus near B*'^B*~ threshold, one might expect 


a{e'^e —)■ B^B^ + X) ~ 1.8 pb , 


( 20 ) 


where X denotes the possibility of one or two additional photons from B* decays. 

The cross section, Eq. (19), may now be multiplied by 2/e(x, ^/s,pr,cut) (see Eq. (2)) 
and integrated over an appropriate range of x. The R* form factor and the 1/s factor in 
the cross section will introduce some suppression, which we shall imitate by introducing a 
maximum Smax = (20 GeV)^. For Eqm = 90 GeV, we thus perform the integral 

a(e'^e~ + X) = — In — / — (35.7 fb) = 1.7 fb . (21) 

71 rrie J X 

Here we have neglected the small deviation of (1 + x^)/{l — x) from 1, taken the limits 
of integration between x^am = (12.72/90)^ and X max = (20/90)^, and used R = 45 GeV. 
One still has to pay the price of the Rc branching fraction to an observable hnal state, but 
as we expect R(Rc —t J/'ipiii') to exceed a percent [65] this seems possible with a sample 
exceeding one event per ab. 

At a different center-of-mass energy, as long as the same range of y/S is taken, one can 
show that the cross section in Eq. (21) scales as (1/s) ln(R/me), so at y/s = 250 GeV, it 
becomes 0.24 fb. 


D Interesting thresholds 

The production of BB pairs has occupied most of the running time of the R factories KEK- 
B and PEP-II. However, some data have been taken at higher energies, as indicated in 
Table H. The GLEO Gollaboration has taken a small amount of data above A;,Ab threshold 
in search of a “magic energy” for Af, pair production; none was found. In Table V we 
summarize some thresholds for heavy flavor production in e’''e“ collisions. 

Here we have used masses tabulated in Ref. [11]. The state R^o in Table V is the 
expected analogue, with = O'*", of the Rso(2317), which is narrow because it lies below 
DK threshold. In order to produce the R^o in e^e“ collisions, it must be accompanied by a 
R* or heavier companion. Angular momentum and parity conservation forbid the process 
^ 7 * —)■ BsqBs. The R^o mass is estimated to be 5717 MeV by assuming that R^o 
and Dsq are chiral partners of Dg and R^ and therefore the R^o ~ Bg splitting is very close 
to the Dso — Dg splitting [66,67]. On the other hand, in order for R^q to be interesting, 
it needs to be narrow. In analogy with DgQ which is narrow because it is below the DK 
threshold, R^o needs to be below BK threshold, i.e., below 5778 MeV. So in any case the 
interesting threshold is between 5717 MeV + rriB* and 5778 MeV + ms*-, i.e., between 
11132 MeV and 11193 MeV. 
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Table V: Some thresholds for heavy flavor production in e'''e collisions. 


Final state 

Threshold 

(MeV) 

BB 

10559 

BB* 

10605 

B*B* 

10650 

BgBs 

10734 

BsB: 

10782 

bib: 

10831 

Bsob: 

11132-11193“ 

AfeAft 

11239 

BcBc 

12551 

Bcb: 

12619-12635^ 

b:b: 

12687-12719^ 

^bc '—'6c 

13842-13890“ 

‘—'66 ^66 

20300-20348“ 


“See text. ^With estimated B*-Bc splitting 68-84 MeV [65]. “ Estimate in [65]. 

VII Conclusions 

While collisions have been studied with impressive statistical power at energies ac¬ 
cessible to the asymmetric B factories KEK-B and PEP-II, the CM energy range from 
about 12 to 80 GeV accessible to PEP, PETRA, and TRISTAN is much less thoroughly 
investigated. The radiative return process —)■ 7 isRe’''e“ —)■ 7isr/, where ISR denotes 

initial-state radiation, can help fill this gap. Some examples are given of processes that 
could be investigated using radiative return, starting from a collider operating at 90 or 250 
GeV. Although the same hnal states / can often be produced with higher cross sections in 
hadronic collisions, the relative cleanliness of the e’''e“ environment gives it an advantage 
whose quantitative value must be investigated using detailed detector simulation. 

Processes which could benefit from radiative return studies at high energies include 
searches for “dark photons” Z', heavy quark (particularly h) production, and spectroscopy 
of states too heavy to be produced at the asymmetric B factories. In studying subenergies 
in the 12-80 GeV range, it was found advantageous to use total GM energies near 
the Z rather than at the highest possible energy. 


Acknowledgements 

We thank Gideon Alexander, Henryk Gzyz, Achim Denig, Stefania Gori, David Tucker- 
Smith, Graziano Venanzoni, and San Pan Wu for helpful discussions. The work of J.L.R. 
was supported in part by the U.S. Department of Energy, Division of High Energy Physics, 
Grant No. DE-FG02-13ER41958, and by funds from the Physics Department of the Uni¬ 
versity of Ghicago. L.-T.W. is supported by DOE grant DE-SG0003930. Monte Garlo 
computations were performed on the Midway cluster supported by the Research Gomput- 


18 







ing Center at the University of Chicago. 


VIII Appendix: Parton luminosity 


Another way to nnderstand the rate for a process below the nominal center-of-mass energy 
is nsing parton Inminosities [68,69]. In this parameterization the cross section is written as 

(j(s) = J dr^^^aabis) , 


where t = s/s and a,b specify the incoming parton species. The hatted qnantities are 
with respect to the colliding partons, primarily electron and photons for colliders and 
qnarks and glnons for pp or pp colliders. The qnantity {1/s){dLat/dr) is called the parton 
luminosity and has nnits of a cross section. It is compnted as 


dLab 

dr 


1 

1 + dab 



fa{x)fb 



+ r (j) A(x) 


(23) 


In the following snbsections we show the parton Inminosity calcnlated for lepton colliders 
and hadron colliders as a means to compare the expected rates at the different types of 
colliders. An important caveat is to note that there are signihcant differences to actnally 
compnting rates at the different machines which means jnst comparing the parton Inminos¬ 
ity valnes can give an inaccnrate pictnre. More realistically one needs to consider leptonic 
versns hadronic branching ratios, detector efhciencies for the hnal states, and backgronnd 
processes. 


A Leptonic parton luminosity 

At an e’''e“ collider when considering the only initial state as e’''e“ and not e’'' 7 , e“ 7 , or 
77 , the electron distribntion fnnction is almost the same as the splitting fnnction fe{x) in 
Eq. (2): 

fe{x,s) = 5{1- x) + -\n — +C>(a^) . (24) 

71 irie \{1 - x)+ 2 J 

The pins distribntion regnlarizes the behavior at a; = 1 [70] and the (3/2)(5(l — x) is for 
overall normalization. The inclnsion of the (5-fnnction is necessary becanse the distribntion 
fnnction is inclnsive and needs to acconnt for the no-splitting case. 

For the parton Inminosity we hnd 


dr 


5(1-r) (l + —\n — + 0{a^)) + —\n 
\ 71 me J 71 


2a , E 1 + 


me 1 — r 


0(cf). 


(25) 


Fig. 10 (left) shows the parton Inminosity for several different CM energies. For narrow 
resonances, we can directly nse the parton Inminosity to compnte rates. According to the 
narrow width approximation for a resonance of spin J, mass m, and total width T 


sd{s) = 47r^(2J -|- l)BiBfmV5{s — m?) , 


(26) 
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Figure 10: Parton luminosities for e+e at y/s = 90 GeV, 250 GeV, 350 GeV (left) and for 
pp at y/s = 1.96 TeV and pp at ^/s = 8 TeV, 13 TeV (right). 


giving a cross section of 

11L \ O LL I / 9/ 

\ / T=m^}, 


(27) 


We can evaluate the appropriate s curve in Fig. 10 at a subenergy m and hnd the rate by 
multiplying by F/m, a spin factor, and branching ratios. As an example, with Fee = 0.322 
keV and m = 10.5794 GeV, we hnd cr(e+e“ —)■ 7 T( 4 S');s = (90 GeV)^) = 9.17 fb, the 
value obtained in Section III. 


B Hadronic parton luminosity 

For comparison we consider a hadron collider with the initial state 

qq = {uu, dd, ss, cc, uu, dd, ss, cc}. (28) 

The parton distribution functions for the proton are non-perturbative functions that de¬ 
scribe the probability a given parton is taken from the proton. 

Fig. 10 (right) shows the parton luminosity for several different GM energies. In con¬ 
trast to lepton colliders, a higher-energy hadron collider also increases the rates at lower 
subenergies. 

The parton distributions are accessed via the LHAPDF interface [71]. The sets used 
are CT10nnlo_as_0118 set [72], MTSW2008nnlo68cl [ 68 ], and NNPDF23_nnlo_as_0118 [73]. 
Only the CTIO set is shown in Fig. 10 as the differences are negligible in the hgure. 

IX Appendix: Fractional luminosity for e+e“ —?► 

In Section IV we discussed fractional luminosity using factorization in the collinear limit 
to compute cr(s) = cr(e+e“ 7 /; s) in terms of a{e^e~ /; s) with the results shown in 
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Figure 11: Fractional luminosity for process —)■ at i/s = 90 GeV with an 

angular cut of 20° (left) and showing the difference between the Monte Carlo and analytic 
calculation (right). A bin size of A = 1 GeV is used. 


Fig. 5. The use of factorization is very convenient to obtain general results independent of 
the hnal state /. 

Here we compute the fractional luminosity exactly with the processes cr(e’''e“ 
and a(e~^e~ to demonstrate that the factorized form used in the main text 

was justihed. We use the Monte Carlo programs MadGraph [58] and Phokhara [74]. Fig. 11 
(left) shows the fractional luminosity compared to our computation from Section IV for 
0^ > 20° and E-y > 2 GeV. The energy cutoff is necessary to cut off the soft divergence. 
Fig. 11 (right) shows the ratio of the Monte Carlo to the analytic factorized form. The 
Phokhara result agrees well with the factorized form and the MadGraph results differ in 
normalization by ~ 5%. Results for other angular cuts are similar. 

Fig. 11 also displays next-to-leading (NLO) order results computed with Phokhara. 
Interpreting these results requires some additional explanation. Recall the dehnition of 
fractional luminosity 

^ L,a{s) , (29) 

CLEcm 

where cr(s) is the three-body cross section and (t(s) is the two-body cross section. 

Explicitly, the left side of Eq. (29) evaluates the differential three-body distribution at 
a given Ecm value and is multiplied by the bin width A to get a three-body cross section 
(i.e., one term of a Riemann sum). The right side of Eq. (29) is the two-body cross section 
evaluated at Ecm — ^ with a coefficient identihed as the fractional luminosity. 

At leading order (LO), this dehnition is unambiguous because the two-body cross section 
is a (5-function in s. At NLO, however, the two-body cross section becomes a distribution 
in s due to real photon emission.il One needs to choose whether to dehne the two-body 

II Considering real emission an iV-body process becomes an -|- 1-body process, but for simplicity we 
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cross section as the integral over all s values or the integral over a small window A near 
the nominal value. 

The former we call the inclusive NLO cross section and evaluate as 



and the latter we call the exclusive NLO cross section 



Both versions are shown in Fig. 11 in green (lower) and orange (upper) dashed lines, 
respectively. 

The inclusive result has a lower fractional luminosity because the two-body cross section 
is larger than the LO result. This is because the real photon emission already induces 
some radiative return. The decrease of the fractional luminosity, relative to LO, at high 
subenergies is due to the three-body cross section decreasing because of the 2 GeV cut on 
photon energy. 

On the other hand, the exclusive result has a lower two-body cross section because the 
integration only includes values of Ecm near the nominal value while the phase space at 
lower values is still populated. The three-body cross section is the same as the inclusive 
case so the same decrease at high subenergies is observed. The exclusive result depends on 
the bin size used. In Fig. 11 we show results using a bin size of A = 1 GeV. As the bin size 
is increased, the exclusive result approaches the inclusive result. 
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